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ABSTRACT 

We present analytic calculations of the electromagnetic torques acting on a magnetic 
neutron star rotating in vacuum, including near-zone torques associated with the in¬ 
ertia of dipole and quadrupole magnetic fields. We incorporate these torques into the 
rotational dynamics of a rigid-body neutron star, and show that the effects of the 
inertial torque can be understood as a modification of the moment of inertia tensor 
of the star. We apply our rotational dynamics equation to the Crab pulsar, including 
intrinsic distortions of the star and various electromagnetic torques, to investigate the 
possibility that the counter-alignment of the magnetic inclination angle, as suggested 
by recent observations, could be explained by pulsar precession. We find that if the ef¬ 
fective principal axis of the pulsar is nearly aligned with cither the magnetic dipole axis 
or the rotation axis, then precession may account for the observed counter-alignment 
over decade timescales. Over the spindown timescale of the pulsar, the magnetic in¬ 
clination angle always decreases. 
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1 INTRODUCTION 

The structure and evolution of magnetic fields is one of the 
key ingredients to understanding various observational man¬ 
ifestations o f radio pulsars and other types of neutron stars 
(NSs ) (e.g., lHarding fe Lail 120061 : iKaspil l20ld : iReiseneggeij 
2013). For radio pulsars, the magnetic inclination angle a, 
defined as the angle between the pulsar’s magnetic dipole 
axis and rotation axis, strongly affects the pulse and po- 
larization profiles in radio and high energy emissions (e.g., 
IRookvard et al.l[2015i ). By analyzing polarization dat a for a 
large number of pulsars, iTauris fe Manchester! (1 19981 ) found 
that, statistically, pulsars with large characteristic ages tend 
to have small magnetic inclination angles, suggesting that 
the magnetic axis align with the spin axis on a time scale 
of order 10 7 year s (see IWeltevrede fe Johnstonl 120081 and 
[Young et al.ll2010l . who found somewhat different alignment 
timescales). On the other hand, general pulsar population 
studies have revealed no evidence for significant torque de¬ 
cay (due to magnetic field decay or alignme nt) over thejml- 
sar l i fetime (~ 10 s yea rs) (e.g., iFaucher-Giguere fe K as oil 

120061: iGullon et al.ll20lll). _ 

Recently. iLvne et alJ (120131 ) found that the radio pulse 
profile of the Crab pulsar has shown a steady increase in 
the separation of the main pulse and interpu lse compo- 
nents at 0.6° per century over 22 years (see also lLvne et all 
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120151) . The increase in pulse seperation was interpreted 
as an increas e in the magnetic inclination angle a. (see 
[Waiters et al.ll2009l ). This interpretation is also consistent 
with departure of the braking index n = (coco)/C j 2 from 3 for 
the Crab pulsar (where uj is the angular rotation frequency). 
Using the braking torque due to a rotating magnetic dipole 
in vacuum, Cj oc — uj 3 sin 2 a, and assuming a constant mag¬ 
netic dipole moment, the braking index is given by 


UJ 

n = 3 + 2- 


a 


uj tan a 


(1) 


With the observed a = 0.6°/century and uj/uj = 

—24. 9 century, and the estimate a ~ 45° dHarding et ali 
mi), we find n ss 2.48, in agree ment with the observed 
value of 71 ~ 2.50 dLvne et al .11 201 .'ll) . 

The increase in the magnetic inclination angle can¬ 
not be explained by the simplest dynamical model of neu¬ 
tron stars (NSs). If one models a NS as a spherical body 
endowed with a frozen-in dipole magnetic field in vac¬ 
uum, one e xpects only a decrease i n the magnetic inclina¬ 
tion angle dDavis fe Goldstein! 1 19701 ). Including the electro¬ 
dynamical effects of the magnetosphere leads to pulsar spin- 
do wn even when the NS has an aligned dipole field (a = 


0) dSpitkovsk j 200d: 


iKalapotharakos et al 


Kalapotharakos_&_^onto 20 ulos 120091 : 
12012 : Tchekhovskov et al. 20131 ) , but 


still predicts magneti c alignment with the rotation axis 
dPhilippov et alJ 12014). While these r esults may b e consis- 
tent with pulsar population statistics dTauris fe Manchester! 


© 0000 RAS 




























































2 J. J. Zanazzi and Dong Lai 


119981) . the short-term (~10 years) increase of a observed in 
the Crab pulsar is unaccounted for. 

A possible physical mechanism for magnetic counter¬ 
alignment (increase of a) is pulsar precession, a topic 
of interest for nearly half a century. Early models of 
free precession modeled the NS as a rigid body un¬ 
dergoing a torq ue due to angular momentum loss from 
dipole radiation (IDavis fe Goldstelnlll97(1 IColclleiclJli97(J'l . 
Then came the inclusion of the pinned superfluid in 
the NS crust, which was sh own to severe l y alter the 
rotat i onal dynamics of NSs dShahaml 19771: Alpar et al.l 
1984l:lAlpar fe Oegelmanlll987l:ICasini fe Montemavorll 19981 : 


Sedrakian et al.1 19991 : Link fe Cutler 20021 ). The effects of 
super-fluidity “destroyed” precession, speeding it up to a 
rate undetectable by observations. But tentative observa¬ 
tional evidence suggested that some pulsars precessed with 
eriods comparab l e to those predicted by free precession 
Suto fe Isol 19851: iTruemper et al.l Il986l: iLvne et all 1 198 A 


IWeisberg et al .IteoioT Makishima et al .1 120141 ). This led many 

to still model the prec e ssion of NSs as free , rathe r than 
forced l| Link fe Epstein| Il997l: iM elatod 1 1 9971. 1 199 91 l200d: 


Jones fe Andersson 200ll : lLink fe Epsteinl 200ll : Wassermanl 


20q 1)7 and to infer interior physics which would lead to weak 

coupling between the crust and the core. Overall, despite the 
uncertainties, free precession remains a possib le mo del for 
understanding the rotational behavior of NSs djone idoU). 

In this paper, we treat the NS as a non-spherical rigid 
body acted upon by electromagnetic (EM) torques. Section[2] 
presents our calculation of the EM torques, including both 
“inertial torques” associated with the inertia of the near¬ 
zone EM field, and dipole radiative torque. In section [3l 
we solve the equations of motion for the NS rotation an¬ 
alytically, reproducin g results which before now were only 
studied numerically dMelatosl 119991 . I2OO0I ). We show that 
the main effects of the inertial torque may be understood 
by modifying the moment of inertia tensor of the NS. In 
section [4] we discuss applications to pulsars, and in partic¬ 
ular to the observed magnetic inclination evolution of the 
Crab pulsar. In section [SJ we summarize our findings, and 
discuss various uncertainties and possible future works. 


2 EM TORQUES ON ROTATING NS IN 
VACUUM 

In this section, we calculate the EM torques on a ro¬ 
tating, magnetized sphere in vacuum. We consider both 
dipole and quadrupole magnetic field topologies. It is well 
known that a rotating magnetic NS must be surrounded 
by a magnetosphere with cu rrent and charge distribu¬ 
tions dGoldreich fe Julianlll969f ). This magnetos phere mod¬ 
ifies the magnetic breaking torque significantly ij Spitkovskvl 
2006; iKalapotharakos fe Contopoulosll2009l : iPhilippov et al.1 
20141 ). However, the torque associated with the near-zone 
magnetic field inertia has not been calculated for magneto¬ 
sphere models. We will show this inertial torque can signif¬ 
icantly affect the precession dynamics of the NS. 


2.1 Dipole Field 

A spherical body with endowed with a dipole field rotating 
in vacuum has two torques acting on it. The first arises from 


the fact that a misaligned spinning dipole emits EM radi¬ 
ation, carrying away angular momentum. We denote this 
torque as T ra( j. The secon d torque arises from the in ertia of 
the dipole magnetic field (IDavis fc Goldsteinll 19701 ). which 
we will denote by Tp. Our calculation yields the expressions 
(see Appendix A) 


_ 2c© , , 

r rad = p-(wxp)xp, 


(2) 


Tp = 


5 Rc 2 


(P'u)(pxu), 


(3) 


where uj is the rotation rate vector, p is the dipole moment, 
p = |p| = BpR 3 / 2, with Bp the magnitude of the fields at 
the magnetic poles, and R is the radius of the NS. 

Note that the numerical coe fficient 3 /5 in f ront of the 
expression f or Tp agrees with [Melatos _Jl997 ), b ut dis ¬ 
agrees with IDavis fe Goldsteinl (Il970l ). Goldreichl dl970l) . 
iGood fc Nel ( 19851) . and~ Beskin et al.l (1 20131 ). all of whom 
quoted slightly different values. Our equation © is ob¬ 
tained by assuming a uniform interior field Bp which ro¬ 
tates rigidly around the spin axis, and an electric field given 
by E = — (v/c) X Bp. Although this interior EM field 

was assumed, equation © only depends on the exterior EM 
field on the surface of the NS [see Eq. (1711) in Appendix]. 

The dif f erence betw een our value and that given by 
IGood fe Nel (Il985l) and iBeskin et al.l (120131 ) may be at¬ 
tributed to the method used to calculate the torque. These 
authors obtained the torque directly through the volume in¬ 
tegral 


T = / r X ( p e E + -j e X B ) dV, 


(4) 


and adopted specific assumptions on the charge den sity p e 
and the current density j e inside the NS. IBeskin et al] (I2OI3I ) 
also included the effects of the EM field’s inertia. The dif¬ 
ference from the value of IDavis fc Goldsteinl (11970i) cannot 
be attributed to such a difference, as the authors appeared 
to h ave used the same me thod to calculate the torque. 
IBeskin fc Zheltoukhovl (120141 ) also obtained the coefficient 
3/5 using the same approach as ours, although they ques¬ 
tioned its validity of the method, suggesting that the final 
answer depends on the internal current distribution inside 
the rotating NS. Indeed, the angular momentum carried by 
the EM field is 

Lem = 4 ^/ rx(ExB)dV. (5) 


The rate of change of Lem due to stellar rotation is of order 


B|Lem| 


BpR 5 fi 2 P 2 n 2 

c 2 Rc 2 ’ 


( 6 ) 


the same order as equation ©• Thus, including Lem into 
the angular momentum equation amounts to a modification 
of the inertial torque expression © by a factor of order 
unity [see Eq. m of section 13.11 . In the remainder of this 
paper, we adopt equation © as the inertial torque and do 
not include Lem into the dynamical equation for the NS. 
Since the precise value of Bp is uncertain by at least a factor 
of two, and more importantly, since the NS suffers a much 
larger deformation than that associated with the field inertia 
(see section [4711) . a correction to equation © by a factor of 
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order unity does not affect the main results of our paper (see 
sections [30. 

For convenience, we define the dimensionless parameter 
tp as 


3 B 2 P R 5 
~ 20 Ic 2 


(7) 


where / is the moment of inertia of NS. With this definition, 
F p becomes 


rp = 7ep(p-o;)(pXu>), (8) 

where p = p/p is the unit vector along the dipole axis. 
Thus, a change of the coefficient in equation m amounts to 
a modification of the value of ep by a factor of order unity. 


2.2 Quadrupole Field 

We have also calculated the inertial torque for NSs with 
an arbitrary magnetic quadrupole moment Q, where Q is 
a symmetric trace-free tensor with eigenvectors qi, q 2 , and 
q 3 , and eigenvalues Q i, Q 2 , and Q 3 (see Appendix). Before 
we state the result, we will explain our field decomposition. 
Because Q is trace-free QT Qi = 0), this quadrupole mo¬ 
ment may be expressed as the sum of two linearly indepen¬ 
dent tensors: Qy and Qi. These two independent trace-free 
quadrupole tensors have components in the {q,} basis given 
by 


(-Q ll/2 

0 

°\ 


Qy = 0 

-Qll/2 

0 1 

(9) 

V 0 

0 

*9ii/ 


IQs 0 

°\ 



Qi = ( 0 — Qs 

0 > 


(10) 


V 0 0 0/ 

where Qy = Q 3 and Qs = (Qi —Q 2 )/ 2 . On the surface of the 
star r = Rr, the radial components of the magnetic fields 
are 


f-B||(f?r) =B||P 2 (f-q 3 ), (11) 

r-B s {Rr) =^ [P 2 (f-q 2 ) - P 2 (r-qi)] , (12) 

where By = —2 Qy/R 4 and Bs = 3 Qs/R 4 , and P 2 denotes 
the Legendre polynomial of order 2. From this, we see that 
the magnetic field of a general quadrupole is completely 
specified by the basis vectors {q;}, and the field strengths 
By and Bs- Figure Q] illustrates the geometry of these field 
components. 

The above composition of the surface quadrupole mag¬ 
netic field completely specifies the external EM fields, which 
give rise to an inertial torque (see Appendix for calculation) 

r Q = Y^-(^-qi)(^xqi) 

4ByBaP J 

+ 525c 2 [(wqaKwXqs) - (wqi)(wxqi)] 

4 R 2 

- 15 4 c 2 [(wqaXwXqs) + (wq 2 )(wxq 2 )]. (13) 

The form of this torque differs from iGood fc Nsl (1 1985h . 
The torque calculated in iGood fc Nel d 1985ll is expressed 
in Cartesian coordinates, with = z, while the expression 



Figure 1. Two components of the quadrupole field. A general 
quadrupole field is specified by three principal axes (qi,q 2 jQ 3 )> 
in which the quadrupole tensor is diagonal, and two elements (Qy 
and Qs) with the corresponding surface magnetic fields J3y and 

B S • 


above is basis independent. For convenience, we define two 
new dimensionless parameters: 


e ll 


e«5 


i b \\ r5 
~ 175 Ic 2 ’ 
4 BfB 5 
= 1575 Ic 2 


(14) 

(15) 


With these definitions, the inertial magnetic quadupole 
torque takes the form 


r Q =/ey(w-qi)(wXq 1 ) 

+2 Iy/e\\es [(u>-q 3 )(o;Xq3) - («-qi)(wXqi)] 

-I eg [(w-q 3 )(o;Xq3) + (w-q 2 )(aJXq 2 )]. (16) 


2.3 Magnetic Inertia Tensor 

To re-express the torques given by equations ([8]) and (1161) , we 
define two tensors Mp and Mg, associated with the inertia 
of the dipole and quadrupole magnetic fields: 

Mp = —I ep(p ® p), (17) 

Mg = -I ey (qr ® qi) 

-21 y/e\\es (q3 ® q3 - qi 0 qi) 

+1 eg (q 3 0 q 3 + q 2 ® q 2 ) . (18) 

The Magnetic Inertia Tensor M is defined to be 

M = Mp + Mg, (19) 

so that the inertial torque takes the form 

Tinert = Tp T Tq = —a; X (M-cl>). (20) 

This new form of the inertial torque Tinert given by equa¬ 
tion ([20)1 will be used in the next section. 


3 NON-DISSIPATIVE PRECESSION 

In this section, we solve the equation of motion for the NS ro¬ 
tation analytically, incorporating a non-spherical rigid body 
moment of inertia tensor and the inertial torques from co¬ 
rotating dipole and quadrupole magnetic fields. We neglect 
the radiative torque, and derive an analytic expression for 


© 0000 RAS, MNRAS 000, 000-000 




















4 J. J. Zanazzi and Dong Lai 


precession period. iMelatod d 19991 . [2000) has previously pre¬ 
sented numerical solutions for such non-dissipative precision, 
incorporating only the dipole torque. 


3.1 Explicit Solution for Non-Dissipative 
Precession 

The Euler equation for the NS rotation takes the form 
duj 

(I-U>) = Tinert, (21) 

where I is the “intrinsic” moment of inertia tensor for the 
NS. The derivative A/At is taken in the frame with the body. 
Substituting expression PHl into the dynamical equation 
m and re-arranging gives 

dcu 

I-— +ux [(I + M).w] =0. (22) 

To solve this system analytically, we define an effective in¬ 
ertia tensor I e ff as 

left = I + M. (23) 

As long as the magnitude |M|/7 is much less than unity, 
or tp, £||, t§ <C 1, we can replace the first term in (1221) by 
Ieff-dw/dt. This approximation is valid in the full regime of 
interest [see Eq. (145II . (14611 . and (1471) below]. Equation (1221) 
then becomes 

^+ti)XL = 0, (24) 

where L = I e fpu> is the effective angular momentum of the 
body, including the inertia term from the near-zone mag¬ 
netic field. 

Equation (1241) is the equation of motion for a freely pre- 
cessing rigid body, and has a well known analytic solution 
which we will summarize. Because I e fi is a real 3x3 sym¬ 
metric tensor, it may be diagonalized. Let e e ff, 3 , e e ff, 2 , and 
e e ff,i denote the three eigenvectors (the principal axes) of 
I e ff with the respective eigenvalues 7 e ff ,3 > l e ff ,2 > leff.i- In 
this frame, equation (1241) has the conserved quantities 

r 2 1-2 r 3 

+ + = 2E , ( 25 ) 

Jeff, 1 -ieff,2 Jeff,3 

L\ + L 2 + L 3 = L 2 , (26) 

where E is the rotational energy of the body and L is the 
angular momentum. The evolution of the components of 
L = L/L can be obta ined by sol ving equation (1241) (see 
lLandau fc Lifshit3l 19691 Chapter 6: lAkgiin et, alJl2006l l: 

L\ = -Acn ((f), k 2 ), 

1/2 = — AvT+l) 2 sn(<(>, A: 2 ), (27) 

L 3 = dn(<f>,k 2 ), 

when L 2 > 2 EI e g t 2 and 

Li = -Adn(k(f),k- 2 ), 

7/2 = — \/(l — A 2 )(l + e -2 ) sn(k(f), k~ 2 ), (28) 

7/3 = y/1 — A 2 cn (kefr, AT 2 ), 

when L 2 < 2 EI e s,2- Here cn, sn, and dn are the Jacobian 


Elliptic functions, and 


, / 7 e ff,i(2757 e ff,3 — T/ 2 ) 

"V 7/ 2 (7 e ff j3 7 eff i) ’ 

(29) 

(p — t UJp 5 

(30) 

CeffT/%/1 — A 2 

7 eff , 3 Vl + e 2 ’ 

(31) 

k 2 - ^ 

_ 1 - A 2 ’ 

(32) 

-Ieff,3 ^ 7 e ff 1 

Ceff = T , 

Jeff, 1 

(33) 

c 2 _ 7 e ff ! 3(7 e ff i 2 — 7 e ff,i) 
7 e ff,i(7 e ff i3 — 7 e ft,2) 

(34) 


Equations (EH imply precession around eej^, and equations 
El imply precession around e e ff,i- This shows the main 
effect of the inertial torque is to modify the equations of 
motion from that of a freely-precessing body with moment of 
inertia I to a freely-precessing body with a modified moment 
of inertia I e ff. 

An effectively biaxial body corresponds to the special 
cases of e = 0 or e = oo. When e = 0, equations H2ZD simplify 
to 

T/i 

U 

7/3 

where 

(j) = COS 9 6 e ff LO t (36) 

with A = sin 9 and cos 9 = u>-e e ff, 3 . When e = oo, equations 
cn simplify to 

Li = — cos 9, 

7/2 = — sin 9 sin(fc<^), (37) 

L 3 = sin# cos(k(f>). 

where 

k(p = cos#e e ffwf (38) 

and A = cos 9 = a>-e e ff,i. 


3.2 Numerical examples 

In this subsection, we illuminate the solutions of equa¬ 
tion 01 with some illustrative examples. We define a 
body axis {ei}, corresponding to the eigenvectors of the 
axis-symmetric tensor I, with an intrinsic ellipticity e = 
(I 3 — Ii)/Ii and a dimensionless angular velocity param¬ 
eter = lj/uj. Figures [ 2 ] and d show the time evolution of 
u>. In both figures, an intrinsically biaxial star is assumed, 
with e = 10 -11 (see sec. m and the symmetry axis along 
£ 3 . 

Figure [2] shows how the dynamics are modified by 
a dipole magnetic field with tp = 2.4 x 10 -12 and di¬ 
rection p oriented in the (13) plane with angle x = 
10°. The fact that e e ff, 3 -c 2 > — const, indicates that cj 
is precessing around e e ff, 3 , with slight variations due to 
the fact that I e fi is slightly triaxial. Figure [3] adds a 
quadrupole magnetic field with eii = tp, es = 0 , and 


= — sin 9 cos <j), 

= — sin 9 sin cj>, (35) 

= cos 9, 
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^3 - ^ef,3 



Time ( tot) 


Ml - CJ e fF,l . ^2 - ^eff,2 



^3 - ^eff,3 



Figure 2. Components of the angular velocity unit vector 
lJ = uj/ uj obtained by integrating equation (1221) . The compo¬ 
nents are Cji = (dashed line) and lj g ff ^ = e e ff pQ) (solid 

line). We assume a biaxial star with ellipticity e = 10 —11 , 
and p oriented in the (13)-plane with angle x — 10° , and 
ep = 2.4 x 10 -12 . We assume that eq = 0. In the {e*} ba¬ 
sis, ei = (1,0,0), e eff) i = (0.9986,0,0.05277), e 3 = (0,0,1) and 
®eff,3 — (—0.05277, 0,0.9986). Note &2 = e e ff,2- The initial condi¬ 
tions are ei*u> = sin(50°), e 2 *w = 0, and e 3 *tj( 0 ) = cos(50°). The 
fact that e e ff 3 *o> ~ const, indicates that Cj is precessing around 
e e ff 3 , and the small variation is due to the fact that I e ff is slightly 
triaxial. 


Figure 3. Same as Figure [2] except we have added a 
quadrupole inertial torque with ey = ep, = 0 , and qi = 
[sin(60°) cos(205°), sin(60°) sin(205°), cos(60°)] in the basis {e^}. 
Because = 0, q 2 and q 3 are irrelevant. Notice that e e ff 3 *a> is 
no longer constant, because the tensor I e ff is much more tri-axial. 


e = 0 or e = oo, these relations simplify to 


Ti =T 3 


2tt 

6 e ftiU COS 0 


(41) 


One can show that the periods (1391) and (1401) give the periods 
displayed in Figures [2] and [3] 


qi = [sin(60°) cos(205°), sin(60°) sin(205°), cos(60°)]. No¬ 
tice that e e ff ? 3 *uj is no longer constant, which follows from a 
more triaxial I e ff than the example displayed in Figure ©■ 
The solutions to equation m are periodic in time, with 
periods Ti and T 3 for precession around e e ff : i and eeff^ re¬ 
spectively. The explicit relation for T\ and T 2 are given by 


Ti = 

4/eff,l 

/1 + e - 2 K(k~ 2 ) 

(39) 

CeffT \ 

1 1-A 2 k ’ 

t 3 = 

4/eff,3 

CeffT \ 


(40) 


where K is the complete elliptic integral of the first kind. 
The reason why we do not have the simple relation T = 
2-k/uj v is because the Jacobi Elliptic functions are not 27r 
periodic, but rather 4 K(k 2 ) periodic. For the special cases 


4 APPLICATION TO PULSARS 

In the previous section, we showed that the effect of the iner¬ 
tial torque is to add an effective deformation in the moment 
of inertia of the NS. In this section we apply this formalism 
to pulsars. We begin with going over the relative magnitudes 
of NS deformations, then examine the behavior of the mag¬ 
netic inclination angle over different timescales, and lastly 
discuss the possibility of explaining the magnetic counter¬ 
alignment of the Crab pulsar with precession. 

4.1 Neutron Star Deformations 

There are several physical causes for the moment of inertia 
tensor I of a NS to depart from spherical symmetry. The first 
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is due to interior magnetic fields, which gives an intrinsic 
body ellipticity t = (I 3 — Ii)/Ii of order 


■ = P 


r 4 b 2 

GAP 


= 2 x 10 -12 ^ 


B, 


10 12 G 


R 


10 6 cm J \ 1 AMq 


M 


(42) 

where B * is the internal magnetic field strength, and /3 is 
a dimensionless constant which depends on the geometry 
of the internal field. For a di pole or toroi dal magn etic field 
topology, p is of order unity IlMastrano et al.l[2013lf . 

The second source is rotation. A uniform density fluid 
star rotating with angular velocity u> has an ellipticity of 

15 c j 2 
167T Gp ’ 

where p is the density. The NS is likely to have deformations 
of this order early in its lifetime. Once the NS crust crystal¬ 
lizes, the body is able to support hydrostatic stresses. If we 
idealize the crust as having a uniform shear modulus p, the 
part of the ellipticity which does n ot align with the rotation 
axis is (iMunk fe MacDonaldlll975h 


Cfluid 


(43) 


^elastic 


1 + p 


Cfluid — 2 X 10 


R 


1 sec. 


10 6 


10 30 dynes/cnT 
7 / \ -3 


M 


(44) 


1.4 M 0/ 

where p = 19p/(2pgR), g is the surface gravity, and p is the 
fiducial value for the shear modulus for the crust, evaluated 
at the density of order ~ 10 14 g/cm 3 . This is a simple order 
of magnitude estimate of e e iastic- More detailed calculations 
for realistic NS parameters may be found in Cutler et al . 
(120031') . 

We compare these “intrinsic” deformations in the mo¬ 
ment of inertial I to the effective moments of inertia induced 
by co-rotating magnetic fields. These are 


tp = 1.5x10 


-13 ( Bp 


V10 12 G/ \ 10 6 cm 


R 


M 


1.4 M q 


(45) 


c II = 5.7x 10“ 


B ,1 


( R 


10 12 G J V. 10 s cm 


M V 
1.4M 0 ) 


(46) 


t s = 2.5x10“ 


B s 


f R 


10 12 G J V. 10 6 cm 


M 


1.4 Mq 


(47) 

Thus, the biggest correction to the moment of inertia tensor 
I e ff comes from tp, assuming that the quadrupole and dipole 
field strengths are similar. In order for the corrections due 
to the quadrupole field to dominate over that of the dipole, 
one needs either By > 5 Bp or Bs > 7 Bp. 

We note that tp is always at least an order of magnitude 
smaller than e mag (since B* > B), and is smaller than e e iastic 
for realistic NS parameters. Of course, tp may be directly 
inferred from the measured P and P of the pulsar: 

(48) 


This gives an observational lower bound on the effective el¬ 
lipticity e e ff associated with I ef r [see Eq. 

l^effl | CP H - ^elastic H" ^mag “h C|| + e<s| > |ep|. (49) 


4.2 Evolution of magnetic inclination angle: 
analytic result for biaxial bodies 


Before presenting general results and applications in secti on 
14.31 we first summarize the key analytic results of lGoldreichl 
(j 1970f ) for an effectively biaxial body (e = 0, 00 ). We define 
three angles y, 9, and a by 


cosy = p-e e fr, 3 , cos# = u>-e e fr, 3 , cos a = th-p. (50) 


We assume that the magnetic field axis p is frozen into the 
body, so that the angle x is constant in time. The other two 
angles in general will evolve. 

On timescales much shorter than the pulsar spindown 
time, the variation of a is due to precession: 


do; 

— « ix!psinxcscasin#sin(u;pt), (51) 

where oj p = cos 9 t e s co. 

Over timescales comparable to the pulsar the spi n down 
time, the precession can be averaged out, giving dGoldreicbl 
[1973): 


1 / did 
uj \ dt 


2 p 2 uj 2 

' 3c 3 / 


sin 2 x + sin 2 9 [ 1 — ^ sin'’ \ 


1 / d sin 9 

sin 9 


dt 


2 p 2 uo 2 
3 c 3 I 


cos 2 Oil-- sin 2 x 


(52) 

(53) 


From relation ([53]). we see that 9 always evolves toward 0° 
or 90°, depending on if the angle \ is greater or less than 
the critical angle Xcrit = sin -1 (-y/2/3) ~ 55°. This evolution 
takes place over the radiative timescale 


3c 3 / 

rrad “ 2^ 

X ( — \ 

X V 1.4 Mq J 



-4 


Bp y 2 

Io^gJ years - 


(54) 

The magnetic inclination angle a is related to 9 and 
(j> = ujpt (the precession phase), and x, through the relation 


cos a = <li-p = cos x cos 9 + sin x sin 9 cos <j>. (55) 

Using equation (1531) and averaging over (j), we obtain 


d sin 2 a 
dt 


2 2 
p lj 

3PI 


sin 2 (2 9 ) 1 - - sin 2 x 


(56) 


This shows that a always decreases over timescales much 
longer than the precession period, but not necessarily to 
zero. From equation l[33|l. we see if 9 —» 0, then a. evolves to 
X, and if 9 —» n/2, then a evolves to n/2 — \- 


4.3 Counter-alignment of the Crab Pulsar 

The Crab pulsar has P = 0.0331 s and P = 4.22 x 10“ 13 s/s, 
implying the characteristic age of P/2P = 1240 years, and 
dipole field of order Bp ~ 4 x 10 12 G. Through modeling 
the Crab pulsar pulse profile, many authors have estimated 


© 0000 RAS, MNRAS 000, 000-000 















































Electromagnetic Torques, Precession and Evolution of Magnetic Inclination of Pulsars 7 


a to be in the ra nge of 45° — 70° (lHarding et al.l 120081 : 
IWatters et al.ll2009l : lDu et al.ll2012l ). The minimum effective 
cllipticity of the NS arises from the inertia of the dipole 
field, tp ~ 2.4 x 10 -12 . This would give a minimum preces¬ 
sion frequency oj v > e v LJ ~ 0.8°/year. The observed da/d t 
is 0.6°/century dLvne et al.ll2013T ). Since da/dt, to ex¬ 
plain the observed da/dt with precession [see Eq. (151 11 1, we 
require either \ 1 (the effective principal axis e e ff ,3 is al¬ 

most aligned with the dipole axis) or 9 <C 1 (befj^ is almost 
aligned with the rotation axis). This would correspond to 
one of two special cases: the NS is dominated by stresses 
from the dipole magnetic field, or from rotation/ elasticity. 

Figures [5] and [5] depict two examples of the evolution of 
the magnetic inclination angle for a NS with an effectively 
biaxial I e g and e e g = 4x 10” 11 . This value could result from 
the elastic part of the rotational distortion [see (1441) 1 or from 
the magnetic distortion associated with an internal field B » 
larger than the dipole field [see m\- Figure [I] shows the 
case with x = 0.15°, so that the principal axis €^3 is nearly 
aligned with p, while Figure [5] corresponds to the case with 
an initial 6 = 0.1°, so that e e ff ,3 is nearly aligned with ah In 
both cases, a increases during half of the precession phase, 
with d consistent with the value 0 . 6°/100 years observed by 
iLvne et ahl (120131 '). 

Figure El depicts an example similar to Fig. [I] except 
with e e ff = 4 x 10 -11 and \ = 1.0°. This value of e e ff is close 
to the lower limit set by tp (associated with the inertia of 
the dipole field). Note that with a mixture of comparable 
inertial poloidal and toroidal fie lds, the magnetic dist ortion 
£ ma g [Eq. (l42l) 1 can be reduced dMastrano et al.ir2013l '). The 
rotational distortion e e i as tic [Eq. (I44D [ would not affect the 
precession if aligned with the spin axis. Both sets of parame¬ 
ters (Fig.[4]and[6]) can account for the observed da/dt of the 
Crab pulsar. In the case of FigureE] the d > 0 lasts for ~ 200 
years because of the long precession period (2n/uj p oc e~f), 
whereas in the case of Figure [4] the d > 0 phase lasts ~ 20 
years. 

In both Figures 2] and El we see a secular decrease of 
a over many precession cycles, with (d) ~ — 1°/100 years. 
Figure El does not display a large secular change in a. This 
is because equations (1521) and E3> predict 6 stays close to 
zero, so a ~ constant according to equation ([561). 


4.4 First and second order braking indexes 

Also plotted on Figure El is the braking index n = ( ujuj)/Cj 2 , 
computed using equation ©■ We see that during the time 
when d ~ 1°/100 years, the currently observed value of 
n = 2.5 is reproduced. However, during the time when 
d < 0, a much larger value of n is expected. Thus continued 
observations of n in the coming decades will test whether 
precession is responsible for the currently observed d for the 
Crab pulsar. 

Modeling the dynamical evolution of the Crab pulsar 
through precession also gives predictions of the second or¬ 
der braking index, defined throu gh m = foiu) 2 )/aj 3 . By re¬ 
writing in terms of n and h rtMelatoall997l ). we have 

m = n(2n — 1) + f-ii, (57) 

OJ 




Figure 4. Evolution of magnetic inclination angle of a NS with an 
effectively biaxial I e ff, with e e ff = 4x 10 —11 . The initial conditions 
are e e ff l'^ = sin(50°), e e f£ 2 ’^ = 0> and e e ff 3 *u> = cos(50°). 
The upper panel shows the variation over decade timescales. We 
assume x = 0.15°, which here is the angle between the dipole axis 
and the principal body axis e e ff 3 . 


and using equations m and ED, 

m = n{2n — 1) + (3 — n)(n — 1) 

— 2(d) 2 ^4-) (cot 2 a + esc 2 a) + 2 ^4-) d up cot (copt). 

(58) 

With the observed d = 0.6°/century and lj/lj = 

—24. 9 century, and the estimate a ss 45° (lHarding et all 
l2008| ). we find 

m x, io.i + 3.5 (yyrr) cot ( w ^)- ( 59 ) 

If cot(oipi) fa 0, we find precession gives an excellent agree¬ 
ment with t he inferred 1993 value of m ~ 10.1 for the 
Crab pulsar ([Lvne et al . Hi). An evolving precession phase 
(u)pt) may in part explain t he increase of the second order 
breaking index to m ~ 45.6 dLvne et alj|2015l) . 

We caution that our model does not include magneto- 
spheric effects on the electroma gnetic tor q ue. F or example, 
if we adapt the spindown law dSpitkovskvll2006l) 



lu oc —tu 3 (l + sin 2 a), 

(60) 

then equation fD should be changed to 



^ n wsma cos a . 

n — 3 + 2--^— a. 

to 1 + sin a 

(61) 

Thus, 

a ~ 

to account for n ~ 2.50 

1.7°/100years, assuming a ~ 45° 

would require 
. The evolution 


of n would be modified, along with the value of m. 
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Time (years) 


Figure 5. Same as Figure [4] except with x — 65° and an initial 
6 = 0 . 1 °, which is the angle between the rotation axis and the 
principal body axis e e ff 3 . In other words, the initial conditions are 
e e s,r<i(0) = sin(0.1°)sin(65°), e effi2 -u(0) = 0, and e eff| 3 -£(0) = 
cos(0.1°) c.os(65°). 




Time (years) 


Figure 6. The evolution of magnetic inclination angle a (upper 
panel) and the braking index n (lower panel). The parameters are 
the same as Fig. [I] except with e e ff = 4 X 10 —12 and \ = 1-0°. 
The currently observed value n = 2.5 is indicated by the dashed 
line. 


5 DISCUSSION AND CONCLUSION 

In this paper, we have studied the rotational dynamics of 
magnetic neutron stars (NSs), modeled as a non-spherical 
rigid body acted upon by electromagnetic (EM) torques. 
First, we presented new calculations of the near-field EM 
torques associated with the magnetic inertia of the dipole 


and quadrupole fields of a NS in vacuum [Eqs. ([3]) and (1131) 1. 
Our analytical results show that if the NS has a quadrupole 
field a few times stronger than the dipole field, the inertial 
quadrupole torque can become more important than the cor¬ 
responding dipole torque. Second, we showed that, despite 
the complexity of the inertial torque expressions, the effects 
of these torques amount to a modification to the intrinsic 
moment of inertia tensor of the star [see Eqs. ©-(usd and 
mi In general, the effective moment of inertia tensor is tri- 
axial even for an intrinsically biaxial NS. This allows us to 
understand analytically the key effects of the inertial torques 
on the precessional dynamics of magnetic NSs. Finally, we 
applied our theoretical results to the Crab pulsar in order 
to understand the physical origin of the recently observed 
co unter-alignment of the p ulsar’s magnetic inclination angle 
q dLvne et al.ll2013l . 120151 ). We showed that it is possible to 
explain the increase of a on decade timescales through pre¬ 
cession. However, since the typical precession frequency co p is 
much greater than the observed rate da/d t, this explanation 
requires some fine-tuning of the principal axes {e e ff,i} of the 
effective moment of inertia tensor I e ff: e e ff : 3 must be almost 
aligned with the magnetic dipole axis or the spin axis. This 
apparent fine-tuning may be expected if the star’s intrinsic 
distortion arises primarily from the dipole field or from the 
rotation. Over timescales comparable to the pulsar’s age, 
the magnetic inclination a always decreases (see Figs. 141(31) . 

If the observed counter-alignment of magnetic inclina¬ 
tion of the Crab pulsar is indeed caused by precession, then 
the sign of da/d t will switch to negative after half a preces¬ 
sion period. Future observations would provide useful test 
and constrain the precession frequency (and thus the distor¬ 
tion of the pulsar). Note that the upper bound on the time 
which the pulsar should take to reverse its counter-alignment 
behavior is t < ~ 200 years [see Eq. (l45l) l. Thus it 

is not unreasonable to suspect that the behavior of a may 
switch from counter-alignment to alignment within a human 
lifetime. 

There are several complications and uncertainties ne¬ 
glected in our model and calculations: 


(i) Our treatment of the crustal elasticity gave only an 
order of magnitude estimate of e e i as ti C , ignoring complica¬ 
tions such as the NS equation of state and the thickness of 
the crust. 

(ii) The internal magnetic field structure was mostly ig¬ 
nored, restricting our evaluation of e mag to an order of mag¬ 
nitude estimate. 

(iii) We assumed that the only dissipative process was 
the radiative spindown to rque. O ther dissipative processes 
such a s polar wande r ing (lMa cy||l974 |) and crust-core cou¬ 
pling s iShuham 197^ AljsaretalJ 19841: Alj3ar_&_Oegehnan] 


19871: Casini fc Montemavoi 19981 : Sedrakian et al.l 19991 : 
Link fe Cutler! |2002| ) may sigificantly affect the NS rota¬ 
tion/precession dynamics. 

(iv) The exterior of the NS was assumed to be a vacuum, 
which is well known to not be true. The presence of a mag¬ 
netosphere may change the inertial torques acting on the 
NS. 


Any one of these effects may affect the dynamical evolution 
of a, and further work is necessary to determine if their 
inclusions will change the main results of this paper. 
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6 APPENDIX: EM TORQUES ON A 
ROTATING MAGNETIC SPHERE IN 
VACUUM 

We begin with the magnetic held in the co-rotating frame 
of the star “frozen in” to the body. In the inertial frame, 
the magnetic held rigidly co-rotates with the star, main¬ 
taining the same shape and magnitudes as in the rotat¬ 
ing frame. If the star is perfectly conducting, the electric 
held inside the sta r in the inertial frame is given exactly by 
E = — (v/c)xB llLandau fe Lifshft3 ll96(J l. If we work in 
the spherical coordinates ( r,9,q i), the coordinate transfor¬ 
mation from the co-rotating frame to the inertial frame is 
cp —» (p — cut, where co is the rotation rate. If the surface of the 
star is given by r = R, the electromagnetic helds must sat¬ 
isfy the boundary conditions (r-B) r =K+ = (r-B) 7 .=H_ and 
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(rxE) r= i{_ = (r xE)jj = jj+. Thus, the external EM fields 
of the rotating star are uniquely determined by the normal 
components of the magnetic field at the stellar surface. 

The complete solution of the EM fields_ i n vacuum in 
terms of multipole moments is derived in Ijacksonl dl998ll . 
These fields, seperated into time-independent terms Bo 
and Eo, and time-dependent terms E' and B', may be re¬ 
expressed in terms of vector spherical harmonics: 


'■£o\ 
x J 


Bo = (£ + 0) ^ 

a M (l, 0) (^) <+2 'E« 3 

0 = X (e + i)a E (e,o)^y +2 
(AQ)p £+2 


- ^0 


En — 


— CLE{ 


x) 


(62) 


(63) 


l,mj^0 


b' = ^2 \ —j 


VW+T) 

+ y/£(£ + l)a M (t,rn) he ( mX ^ ' 


+ 


cim{£, 


he(mx) + mxh' e (mx) 


y/l(l + 1) L mx 

e' = ^2 | - vw + i ) a E (a m ) hl 




(64) 


t lm 


cle(£, m) 

y/l{l+l) L 

. Qm(£, m) 


mx 

ht{mx) + mxh' e (mx) 


mx 




y/W+rj 


ht(mx)<&lm 


where 


Enn — y£m(@ , kjt)f, 

^>em = rVYe m (6, <j> - ojt), 

= r xVYi m (6,tf> — ut), 


(65) 


( 66 ) 


with x = ruj/c and Y) m denote spherical harmonics. The 
actual fields are given by the real parts of E and B. Defining 
qi m = J dfl(B ■ r) r =R, the boundary conditions give the 
multi-pole moments: 


a M (£,m) 


a E (l, m) 


I i +i — u 

e(t+i) [(1 — $to)£J(e+i)oq(e+i)o 
— (£ + 1)-/«o9(£-i)o] 

xq _ mxQ _ 

< yJ7(l+ 1) h( - ('mx 0 )+mx 0 h' £ (mx 0 ) 

x [£J(e+ 1 ) m Q(t+ 1 ) m 

m] 


(67) 

m = 0 


m yf 0 

( 68 ) 


where xo = Ruj/c , Stm is the Kronecker delta, and 


Jim — 


l 2 —m 2 
Al z — 1 

o 


if |m| < ■ 
if Iml > 


(69) 


The Kronecker delta term forces charge neutrality. 

It may be shown th at for l = 1 the abov e expressions re - 
produce the solution bv lDeutschl (Il955fl [see iMelatosI (Il997l ~) 


for corrections]. For this work, we are interested in both the 
t = 1 and t = 2 terms. Let the unit vector in the direction 
of the dipole moment be p, and the three unit eigenvectors 
for the quadrupole moment be qi, < 52 , and q 3 (see Fig [I]). 
If the magnetic field consists only of £ = 1 term, we may 
define the magnitude of Bp as Bp = B(R, p), where the 
notation B(r,h) denotes the magnetic field evaluated at r 
with angular coordinates in the direction of the unit vec¬ 
tor n. The quadrupole magnetic field has two independent 
components, due to the fact that the quadrupole tensor Q is 
trace free (see Fig.[l]). We define these to be By = B(R, qi) 
and B s = \B(R,q 3 ) -B(B, q 2 )]/2. 

The time-averaged EM Maxwell stress tensor is given 
by 



E <g> E* + B <g> B* 


i<5(|E| 2 + |B| 2 ) 


(70) 


where 8 is the unit dyadic, and the stars denote the complex 
conjugate. One may decompose T into spherical coordinates 
to evaluate the torque T in Cartesian coordinates: 


Ta, = —R J dO [T r g sin <j> + T rr f, cos 0 cos 4>\ r=R , 

= R f dff [T r g cos 4> — T rt t, cos 6 sin <j>\ „ , (71) 

T z = B/d9[T^sine] r=fl , 


where x + iy = r sin de 1 ^, z = r cos 6. 

After integrating equation CD, we expand the expres¬ 
sions in xo and only keep terms proportional to a: 2 . We then 
arrive at equations © and m- 
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